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Electronic Structure of Stripes in Two-Dimensional Hubbard Model 

Masanori IchiokaQ and Kazushige MachidaE*! 

Department of Physics, Okayama University, Okayama 700-8530 
(Received June 7, 1999 ) 

Focusing on La2-a;Sra:Cu04, we study the stripe structure by the self-consistent mean-field 
theory of the Hubbard model. By introducing the realistic Fermi surface topology, the SDW- 
gapped insulator is changed to metallic. The solitonic features of the stripe structure and the 
contribution of the mid-gap states are presented. We consider the band dispersion, the local 
density of states, the spectral weight, and the optical conductivity, associated with the solitonic 
structure. These results may provide important information for the experimental research of 
the stripe structure, such as the angle-resolved photoemission experiments. The "Fermi surface" 
shape is changed depending on the ratio of the incommensurability 5 and the hole density rih. 
In real space, only the stripe region is metallic when (5/nh is large. 
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§1. Introduction 

The concept of the stripe structure is quite versatile 
and powerful in understanding various spatially modu- 
lated orderings ranging from insulator to metal, such as 
spin or charge density waves (SDW or CDW) or spin 
Peierls systems. A canonical example of such a stripe 
description is the incommensurate SDW in metallic Cr 
dilutely doped with Mn or V atoms. Various aspects 
of the SDW under doping are coherently describable in 
terms of the stripe.tiB Another canonical example for 
such description is the spin Peierls systems such as in 
CuGeOs or other organic materials under a magnetic 
field, wluch effectively acts as changing the electron fill- 
ingE'QB' Thus, it is believed that the stripe concept is 
a useful paradigm for low dimensional electron systems 
with nearly half- fillings. 

Anomalous electronic properties in under-doped high 
Tc cuprates attract much attention because understand- 
ing of the unusual ground state may lead us to a clue 
to a high Tc mechanism. One of the key features is 
the magnetic incommensurate (IC) structure, i.e., the 
stripe structure. Recently, the static magnetic IC struc- 
ture is observed by the -elastic neutron scattering on 
La2-^,;Sr^Cu04 (LSCO).!'^'!) It suggests the static 
stripe structure. The static IC structure is also reported 
on Lai.6-2;Ndo.4Sr2;Cu04 [ref. |l^ and La2-xSr^Ni04+j, 
[refs. |l| and ^fi) In YBazCuaOr-i, the IC fluctuations 
are reported by the inelastic neutron expejira£ints.E3'Ej^ 

Almost a decade ago, one of the authorsEj'Ej^ predicts 
some of the features of the above static IC spin mod- 
ulation in lightly doped cuprate systems, stressing the 
charged stripe or the solitonic structure as a convenient 
and universal "vehicle" to accommodate excess carriers 
in an otherwise commensurate antiferromagnetic (AF) 
state within a mean-field treatment for a simple two- 
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dimensional Hubbard model. Some features q£ 

indenpridputly found by Others at that timeO 

later .tJ^l2i|) And they are now confirmed by a more sophis- 
ticated method, such as DMRG.cS' But, since the above 
solitonic structure is an insulator in nature while high Tc 
cuprates are metallic, it was claimed that the other sce- 
nario such as frustrating charge segregation is needed. Ej-* 
However, we can consider the metallic structure in the 
same framework of the selfconsistent mean-field theory, if 
the realistic Fermi surface topology is taken into account 
by introducing a next nearest neighbor site hopping. It 
is still important to study the detailed structure of the 
stripe state in this framework. The stripe state in the d- 
p model (i.e., including 0-site contribution) was studied 
within the mean-field theory in ref. 

The purpose of this paper is to analyze the stripe struc- 
ture based on the mean field method of the Hubbard 
model. Needless to say, the mean field theory which sat- 
isfies variational principle for the energy minimum can be 
justifiable when a long-range order exists, such as in the 
present case. Thermal or quantum fiuctuations can be 
treated perturbatively from this solution, if it is needed. 
We consider the doping dependence and the effect of the 
realistic Fermi surface topology. Our numerical method 
is described in The profile of the stripe structure and 
its band structure are studied in We also discuss the 
doping dependence of the incommensurability. In we 
study the density of states (DOS) and the local density 
of states (LDOS). The latter will be observed by scan- 
ning tunneling microscopy (STM) experiments. In 
we discuss the spectral weight, i.e., fc-resolved DOS. Our 
results are compared with the data of thrvjjiiagilejresolved 
photoemission (ARPES) experiments,E3E3'E3'E§' which 
is a powerful method to ptiidy Fermi surface topology 
and pseudo-gap structure.E^ In we calculate the op- 
tical conductivity, which detects excitaijians-of the stripe 
structure, and compare with the data.E2l'E2|) Our theo- 
retical study shows that these experimental methods can 
give us a lot of useful information of the stripe structure. 
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such as the mid gap state or the fe-dependent SDW gap. 
The last section is devoted to summary and discussions. 
Throughout this paper, we set h ~ k-Q — \ and the lat- 
tice constant c = 1. Some of our results were briefly 
reported in ref. ^ 

§2. Selfconsistent Mean-Field Theory 

We start out with the standard Hubbard model on a 
two-dimensional square lattice and introduce the mean 
field 

{n,^a)^\{n^^oM^) (2.1) 

at i-site, where (Af^) is the charge (spin) density and 
i = {ix,iy)- Thus, the one-body Hamiltonian is given 

by 

^ = - hjal^Qj^^ 

+C/^(n,,_,K, - C/^(n,,;)(n,,-r). (2.2) 

where Ui^a = al^^ai^a (cr = 1 for t, -1 for i) and 
(oi^o-) is a creation (annihilation) operator. For the near- 
est neighbor (NN) pairs tij = t. For the next NN 
pairs situated on a diagonal position in a square lattice, 
tij = t' . For the third neighbor pairs, which are situated 
along the NN bond direction, tij = t" . For the stripe 
structure with the N site periodicity, 

M= c"«"'(";Q,-)' (2-3) 

0<i<JV 

where is the position of i-site. The modulation vector 
is given by Q = — 2tt{^ — S,-^ — (5) for diagonal in- 
commensurate (DIG) state, Q = Qy = 27r(i, ^ — S) for 
vertical incommensurate (VIC) state, and Q = Qq = 
27r(i, i) for commensurate antiferromagnetic (C-AF) 
state. We set the incommensurability 5 ~ M/N (ratio- 
nal fraction) with even integer A^. We consider cq. (2.2) 
in the reciprocal space. There, k is coupled to k + IQ 
in the term UY,kj,A'^iQ-'y)0'k+iQM°'k.<y "^here a^^^ is 
the Fourier component of ,j. We write k = k^ -\- mQ 
{m — 0,1, ■ ■ ■ , N — I), where feg is within the reduced 
Brillouin zone of the size {2t:)'^/N. Then, eq. (2.2) is 
reduced to 

n = ni + n2 (2.4) 

0</<A' 

= ^fco,<T,a7l„,<,,„7fco,'T,a, (2.5) 

ko ,a.rn 

0<i<JV 



where = (nj^o.o-) is the average number of a-spin 
electron per site, Nk = I]i 1 = J2k„.rn 1' and e{k) = 
eo(fc) + ei(fc) with eo{k) = —2t{coskx + cosky) and 
ei(fc) = — 4t' cos fcx cos fcy — 2t"{cos2kx + cos2ky). In 
eq. (2.5), we diagonalize N x N matrix Hi by using the 
unitary transformation 

(2.7) 

a 

where a(= 1, 2, • • • , N) is the label of the eigen-energy 
^ko,(7,a and the wave function Wfco^o-,Q,m- 
The self-consistent condition is given by 

{niQ,a) = Yl (^ko+rnQ,a'^ko + {m+l)Q,a) 

kQ,m 

^ X] '^ko,cr,a,ni'^ko,cr,a,ni-\-lf{Eko,a,a)- 

kQ,m,a 

(2.8) 

Here, we use (7l„,^,„7fco,'T,a') = Sa,a' f{Ek„,a,a) with the 
Fermi distribution function f{E) = (e^^"'^)/^ + 1)^1. 

In our calculation, the eigen-energy and the wave func- 
tions are obtained by eq. (2.5) under the given (?1;q.(j). 
The chemical potential fi is determined by the condi- 
tion h-^ — — ^{1 ~ nil) for 1 — rih filling. Then, we 
calculate {niQ u) by eq. (2.8), and use it in the next 
step calculation of eq. (2.5). This iteration procedure 
is repeated until a sufficient self-consistent solution is 
obtained. Among the possible stripe structures (C-AF, 
Die and VIC with various incommensurability 5), the 
stable structure is determined by the total energy {{H)) 
minimum at zero temperature. We mainly consider the 
case t' = —^t and t" = 0. This set of the parameters 
is proposed for LSCO to reprodiiqe the Fermi surface 
obtained by the band calculation.E3'EJ'E3 

§3. Spatial Profile and Band Structure 

For t' = Q, since the Fermi surface is almost a square 
shape near half- filling (nh ~ 0), the complete nesting oc- 
cura^and one dimensional (ID) picture may be applica- 
ble.H By introducing the effect of t' , the Fermi surface 
deviates from the square shape. At first, we consider 
the i'-effect on the C-AF state. The energy dispersion 
is split into the upper conduction band and the lower 
valence band, and it is given by 

ii;±(fc) = i[e(fe)+6(fc + Qc) 

±{(e(fc)-6(fc + Qc))' + 4|AAF|'}^/'] 
= ei(fc)±{eo(fc)' + |AAFn'/', (3.1) 

with the AF gap Aaf. The AF gap opens at feg — 
(fca;,±7r — kx) and (/ca:,±7r -I- kx) by the nesting of Qq. 
Then, the lower and upper edges of the AF gap are 
written as E±{kQ) = ei(feQ) ± Aaf. For t' < 0, 
E±{kQ) is small near (1,0) and (0, 1) in reciprocal lat- 
tice units compared with near {^,^), while the gap 
E^lkq) — E^{kQ) — 2Aaf does not depend on the k- 
point. This dispersive AF gap structure is in fact ob- 
served by ARPES on some parent compounds of High 
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Fig. 1. Dispersion relations for DIG along the perimeter of the 
reduced Brillouin zone for njj = 1/16 = 0.0625, U = 3.6t, <5/nii = 
1/2. t' /t = (a) and —1/6 (b). The Fermi energy Ep is at zero. 
Energy is scaled by t. 



Tc supercondiLctors such as a half-filled Mott insulator 
CazCuOzClz-B) 

Next, we consider the i'-effect on the IC state for fi- 
nite nil- We also find the dispersive gap structure in the 
SDW gap. We show the energy dispersion within the 
reduced Brillouin zone in Fig. 1 for DIG and in Fig. 2 
for VIC in the case rih = 1/16 = 0.0625. There, a large 
SDW (or AF) gap opens between the conduction band 
and valence band, and there are some mid-gap bands 
in the large SDW gap. In the figure, unit vectors of 
the reduced Brillouin zone are given by Ui = (0, 27r), 
U2 = {2tt/N, 2tt/N) for DIG, and Ui = (tt, tt - 2nM/N), 
U2 — {0,4:'k/N) for VIG. The horizontal axes in Figs. 1 
and 2 are along the path + M2) - - o 

+ 1X2). 



In the stripe structure with N/M periodicity (i.e., M 
periods in N sites), the size of the Brillouin zone is re- 
duced to 1/A^, and the energy dispersion is split to N 
bands as seen in Figs. 1 and 2. If it were half- filling, 
lower N/2 bands would be occupied by electrons. In the 
case of finite doping rih = M/{N/2), the induced holes 
enter M bands among N/2 bands. Then, an SDW gap 
opens above the fully occupied N/2 — AI bands, and 2M 
bands above the filled bands form a mid-gap state. Since 
the gapped state minimizes the total energy, the stripe 
structure with N/M periodicity is stable in an insulator 
state. Then, the incommensurability is exactly given by 
S = M/N = irih in this insulator state. This relation 
between the incommensurability 5 and the hole density 
nh is confirmed by our energy estimate at T = in the 
insulator state, even for finite t' . This is a generic state- 
ment which might be iiapjQctajit^to .interpret the neutron 
diffraction experimentsMQS&Byy) 

For t' — 0, the lower and upper edges of the SDW gap 
are almost flat [Figs. 1(a) and 2(a)]. But, for finite t' , the 
gap edges become dispersive [Figs. 1(b) and 2(b)-(d)]. 
With increasing \t'\, both the valence band and the mid- 
gap state touch the Fermi energy Ep in the VIG case 
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Fig. 2. Dispersion relations for VIC along the perimeter of the 
reduced Brillouin zone for rih = 1/16 = 0.0625, U = 3.6t. (a) 
t' /t = and (5/nh = 1/2. (b) t' /t = -1/6 and S/uy, = 1/2. (c) 
t'/t = -1/6 and <5/rah = 3/4. (d) t' /t = -1/6 and <5/nh = 1. 
The Fermi energy E-p is at zero. Energy is scaled by t. 



with S = ^rih- And the stripe state becomes metallic 
[Fig. 2(b)]. We set = in the figures. For smaller U, 
it becomes easily metallic since the SDW gap is small. 

The total energy of DIG is_lpwer than that of VIC 
for higher U and smaller nh-tSHEj^ For example, the 
DIG is stable for rih < ric ^ 0.1 in the case U = 3.6t 
and t' = — With raising (lowering) U, ric increases 
(decreases). When \t'\ is larger, Uc decreases. The DIG 
state is difficult to become metallic since the dispersion 
of the mid-gap state is flat as shown in Fig. 1. While 
the metallic DIG state can be realized for larger rih and 
smaller U , the VIC state is stable in these parameter 
regions. In the parameter region of the stable DIG, it 
is an insulator. Therefore, it is unlikely that the DIG is 
realized to be metallic upon doping. This is a generic 
statement within the present framework. 

On the other hand, for the electron doping case (nh < 
0), the G-AF state has lower energy than those of VIC 
and DIG for small jnhj and lower U. And the VIC is 
stable for large jnhj and higher U . The DIG state does 
not appear for nh < 0. The stable VIC's parameter 
region becomes narrower with increasing \t'\. The VIC 
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for nh < is an insulator, because the top of the mid-gap 
band is flat as seen in Fig. 2(b). E-p is located in the 
gap between the mid-gap band and the upper conduction 
band. 

Once the metallic state is realized for nh > 0, the 
relation S = ^rih, which is valid for the insulating case 
due to the energy gain by the SDW gap, is not obeyed. 
Figure 3 shows the Uh dependence of 6 for VIC by our 
estimate of the total energy minimum. There, t' /t is 
changed from to —0.2, and t" — 0. It is found that 6 is 
fixed to be in insulator state {t'/t = 0, —0.05). But, 
in metallic state of VIC, 6 becomes larger than the ^n^- 
line, and approach the rih-line with increasing \t'\. The 
nh-dependence is almost linear at small rih, and slightly 
suppressed at large nh. The nh dependence of S dose not 
depend on U in the metallic state for 3 < J7/t < 4. It 
only depends on t', i.e., the shape of the Fermi surface. 
It may be related to the nesting of Q vector. Figure 
4 shows the Fermi surface of e(fe) for nh = 0.125 and 
the Qv'Shifted Fermi surface. The nesting seems to be 
better in the case 6 ~ |nh than in the case S ~ ^nh. It 
is the reason why 6 is larger than ^nh- With increasing 
nh further, the parameter S/rih for the best nesting is 
shifted to smaller. Then, the saturation of S occurs at 
larger nh as seen from Fig. 3. We also consider the 
case t' — — 0.12i and t" = O.OSt. These parameters are 
proposed to reproduGe-jbhejEermi surface obtained by the 
ARPES experimentsJ23'E§E3' The obtained S is near the 
nh-line, as shown by filled circles in Fig. 3. We note that 
the C-AF state has lower energy than VIC and DIC for 
small nh in the case U is small and (5/nh approaches 1. 
As a result, 5 = ^nh for insulator DIC and VIC. And for 
metallic VIC, 5 approaches 5 = nh with increasing \t'\ or 
\t"\. As for the DIC case, 6 = ^n^ in the insulator state. 
In the metallic DIC, which is difficult to be stable, S is 
suppressed than inh-hne. 

Some of our results are consistent with the experimen- 
tal results. The elastic and inelastic neutron scattering 
experiments suggest thai the insulator La2-a;Srj;Ni04+j, 
is DIC with 6 ~ inh-0'tJ't3^ The metallic (or supercon- 
ducting) iSCO Jor X > 0.06 and Lai g_2.Ndo.4Sr2.CuO4 
are VICfl'LlB&El' There, it is suggested that 5 ~ nh 
for small nh- The saturation of S occurs for large nh, 
where the inter-stripe distance is close to the width of 
each stripe as shown later in Fig. 5(b). The upper 
limit of S seems to be determined by the width of each 
stripe, beyond which S saturates. The LSCO exhibits 
the transition from metallic to insulator at a; ^ 0.05. 
With this transition, VIC with ^ ~ a; (~ nh) is changed 
to DIC with S ^ [in the notation of ref. ^, Qq = 
27r(i — f J ^ ~ f ) with e ~ nh]. Reference ^ reported that 
e = 0.06 ± 0.005 at a; = 0.05 in DIC. Also in our model, 
VIC is changed to DIC at critical nh with decreasing nh 
for appropriate U. The DIC is insulator with 6 = ^nh, 
and VIC can be metallic. Our model shows that 6 be- 
comes larger than ^nh in the metallic VIC case. But, 
it fails to give the relation exactly 5 = nh- To improve 
this point, we need to carefully adjust the parameters of 
the Fermi surface topology.EZP On this problem, the ef- 
fect of the long range Coulomb interaction was also sug- 



0.1 




0.0 



0.1 



Fig. 3. The dependence of the incommensurability 5 to mini- 
mize the total energy for VIC. Open symbols with solid lines are 
for t'/t = 0, -0.05 (insulator), -0.1, -1/6. -0.2 (metal). Filled 
circles with a dotted line are for t' = — 0.12t and t" = O.OSt. 
The lines are guide for the eye. In metallic case, 5 is enhanced 
with increasing \t'\. We also show lines for 5 = nh and 5 = nh/2 
(dashed lines). 
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Fig. 4. Nesting features: the original Fermi surface of e(fc) [solid 
lines] and Qy-shifted ones for 5/n^=l/2 [dotted lines], 3/4 
[dashed lines] and 1 [dash-dotted lines], n^ = 0.125, t'/t = —1/6, 
t" = 0. The 27r X 2-k Brillouin zone is presented. The nesting is 
seen to be better for <5/nh ~ 3/4 than 1/2 or 1. 



gestedra' There is a possibility that the VIC state with 
(5 = nh becomes more stable by adding the long range 
Coulomb interaction to the Fermi surface topology effect. 

To sec the effect of the 5 shift on the gap structure at 
fixed nh, we show the energy dispersion for 5 = |nh (c) 
[which gives the total energy minimum] and 5 = n\i (d) in 
Fig. 2. They show qualitatively the same band structure 
as in the case 5 — ^nh [Fig. 2(b)]. But, the position of 
E-p is changed depending on 5. With increasing 5 at fixed 
nh (i.e., effectively decreasing nh compared with 5), Ep 
is shifted toward higher energy in the band structure. 
Both the valence band and the mid-gap state cross Ep 
for S = inh. For S = |nh and 6 = nh, only the mid-gap 
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(b) ',-site 

Fig. 5. The spatial profile of the charge density rii and the mag- 
netization {—iy^^''ymi along the j/-direction in VIC. U = 3.6t, 
t'/t = -1/6. (a) rih = 1/48 = 0.021, 6/rn, = 3/4. (b) 
nh = 1/8 = 0.125, <5/"h = 5/7. 



state crosses Ep, and the valence band does not appear 
at E-p . In the following, we consider the optimized S case 
shown in Fig. 3 ((5 '--^ |nh for t' — —-^t). But the result 
is not qualitatively changed for other 6. To know other 
S cases, it is enough to shift Ep in the following results, 
keeping the overall gap structure unchanged. 

Next, we discuss the spatial profile of the stripe struc- 
ture. In Fig. 5, we show the spatial profiles of the 
charge density rii and the magnetization (— 1 
along the y-direction. The excess carriers are located 
in the stripe region, and form CDW. The width of the 
charged stripe is about 3 sites in this case. The AF 
changes its phase by tt when it crosses the stripe region. 
The envelope of looks like sinusoidal wave for large rih, 
but it looks like square wave for small n^. To clearly see 
this profile change, we consider the Fourier components 
riiQ = niQ^^ + niQ^^ and rn,Q = ^ - njg |. In Fig. 
6, we plot the rih-dependence of rriiQ/mQ {I = 3,5,7) 
and niq/mq {I = 2,4,6). For large rih, rriiQ/mQ 
(/ > 3). It means that the fundamental component 
ruQ is dominant and that nii varies as a simple sinu- 
soidal wave. On the other hand, with decreasing rih, 
miQ/mq (-l)^'-!)/^ > 3: odd integer). It is the 
Fourier component of the square wave. The continuous 
change from the sinusoidal wave to the square wave with 
decreasing nh is reminiscent of the Jacobi elliptic func- 
tion sn{x,k), which was thcj_£orm in ID CDW or SDW 
systems near the half-filling O 

This TT-h-dependence of the m^-profile is expected to 
be observed by several experiments. We show the n^- 
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Fig. 6. The njj-dependence of the higher harmonics of the Fourier 
components, tti^q (a) and rijQ (b) for VIC (o) and DIC (o). We 
plot ^iqI^q (l = 3,5,7) and ^iq/'i^q {I = 2,4,6), respec- 
tively. Lines are guide for the eye, and • shows the extrapolation 
for rtij 0. U = 3.6t, t'/t = —1/6. We use the optimized value 
of 5 given in Fig. 3 
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Fig. 7. The nij-dependence of the maximum (max|mi|) and the 
variance {{rn^)^^^) of the magnetization in VIC (o) and DIC (o). 
U = 3M, t'/t = -1/6. 



dependence of the maximum max|TOi|, and the variance 
(m^)^/^. As shown in Fig. 7, max|mi| sustains a con- 
stant for < rih < 0.06, but {m^)^/'^ decreases gradually 
with increasing rih- That is, in the experiments which 
detect the maximum of the magnetic moment such as 
the neutron scattering, the magnetization does not de- 
pend on Tih for small rih- But, in the experiments which 
detect the spatial average of the magnetic moment such 
as NMR, the magnetization gradually decreases wtt^ i^" 
creasing nh. This distinction is observed in Cr.tlB As 
seen in Fig. 8, the distribution function P{m) of the 
magnetic moment (i.e., the resonance line shape of NMR 
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Fig. 8. The distribution function of the magnetic moment, P(m), 
for rih = 1/48 = 0.021 {•) and for rth = 1/8 = 0.125 (o) in VIC, 
which correspond to the cases of Figs. 5(a) and 5(b), respec- 
tively. Lines are guides for the eye. 
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Fig. 9. The density of states N{E). U = 3.6t, t' /t = -1/6, 
nh = 1/16 = 0.0625. (a) The insulator DIG case of Fig. 1(b). 
S/ni, = 1/2. (b) The metallic VIC case of Fig. 2(c). S/ni, = 3/4. 
The Fermi energy Bp is at 0. 



or /uSR) changes its shape reflecting the change from the 
square wave to the sinusoidal wave. At smaller rih, P{m) 
splits into two sharp peaks around ±max|mi| as in the 
AF state. At large rih, P{m) distributes broadly between 
the two peaks. 

As for the profile of the charge density, in the limit 
of large rih, n2Q/mQ is finite and other n;Q — > > 
4). It means that is reduced to the simple sinusoidal 
wave with half period of m^. In the limit of small rih, 
riiQ ~ — (— l)'/^nh ^ (/ > 2: even integer). It reflects 
that the excess carriers are restricted in the (5-function- 
like stripe region whose width is a few sites, and that the 
other wide region is the insulating half-filled AF state. 

The spatial profile of the stripe structure obtained here 
is qualitatpiply, |thc same as that of the insulator state 
at t' = Omi3H> Our results show that the solitonic 
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Fig. 10. Local density of states N(r, E) in the metallic VIC case 
of Fig. 9(b). (a) At the center site of the stripe, (b) At the 
nearest neighbor site to the center site, (c) At the next nearest 
neighbor site to the center site, (d) At the farthest site from 
the stripe. With leaving the stripe [(a) to (d)], the contribution 
from the mid-gap state disappears. Dashed lines show the up- 
spin contribution. U = 3.6i, i'/i = -1/6, rih = 1/16 = 0.0625, 
'^/'^h = 3/4. The Fermi energy E-p is at 0. 



profile in the insulator state remains unchanged even in 
the metallic state. But the electronic state is changed in 
the metallic state. We discuss it in the following sections. 



§4. Local Density of States 

In terms of the eigen-energy and the wave-function, 
Green's function can be written as 

G.iry.E) ^ "'r:"i''^"l°'"°^''^ (4.1) 



where 



= N;"'Y.e<'^-+"^Qy-^uu,,^^^^^ (4.2) 



and rj{> 0) is an infinitesimally small constant. We ob- 
tain the LDOS by 

N„{r,,E) = --lmG„{r,,r,,E) 

TT 
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Fig. 11. Contour map of the spectral weight N{k, E) along the symmetry lines in the original Brillouin zone for the metallic VIC case 
of Fig. 2(c). (a) Along the line (0, 0) - (1, 0) - (1, 1) - (0, 0) in reciprocal lattice units, (b) Along the line (0, 0) - (0, 1) - (-1, 1) - (0, 0). 
(c) Along the line (-1,0) - (0, 1) - (1,0). U = 3.6t, t'/t = -1/6, rih = 1/16 = 0.0625, 5/nh = 3/4. The Fermi energy is at 0. 
Energy E is scaled by t. 



= E Wko,.Ari)\^^{E-Ek,,,^^). (4.3) 

ko,a 

The DOS is the spatial average of LDOS, i.e., 

i ko,a 

(4.4) 

Figure 9 shows the DOS for DIG and VIC in the case 
rih = 1/16 = 0.0625. The mid-gap state is seen within 
the large AF gap. In the insulator DIG case, there is 
a small gap between the mid-gap states and the valence 
band, and is located in the small gap. In the metallic 
VIG case, the small gap is buried in this DOS plot. For 
the case 5 = ^nh, E-p{= 0) is located at the overlapped 
mid-gap state and the valence band. With increasing 
6 at fixed rih, E-p is shifted to higher energy without 
qualitative change of the DOS structure. But, in the 
large (5(~ rzh) case for large nh(> 0.1), the mid-gap state 
and the valence band are separated in the DOS plot, 
where Ey is located in the mid-gap states. As a result, 
even in the metallic state, a gap feature in the electron 
state still remains. This point is discussed further in the 
next section. 

Next, we consider the LDOS, N{ri,E). It is shown 
in Fig. 10. The contribution of the mid-gap state is en- 
hanced at the center of the stripe, and decreased with 
leaving the center. Far from the stripe, the LDOS is like 
an AF-gapped spectrum since the mid-gap state's contri- 
bution vanishes. For 5 = |nh (and also for (5 = nh), E-p 
located in the mid-gap state is in the AF gap outside the 
stripe region [Fig. 10(d)]. Then only the stripe region 
is metallic, and the other AF region is an insulator. If 
(5 = ^nii in the metallic VIG case, where Ep is effectively 



shifted to lower energy, Ep touches the mid-gap state in 
the stripe region and the valence band in the outside 
region. Then, both the stripe region and outer AF re- 
gion are metallic. If these characteristic features of the 
r-resolved DOS structure arc observed by experiments 
such as STM, it may be a powerful means to confirm the 
solitonic features of the stripe structure. 

§5. Spectral Weight 

In this section, we study the fe-resolved DOS, i.e., spec- 
tral weight at each fe point. It is given by 

N^{k,E) = --lmG^{k,E) 

TT 

= E \^ko,a,a,mfHE - Ek^^„^JS{ko + mQ - k), 

(5.1) 

where 

G,{k, E) = E e-'*"('-'-'^C;<.(r, r', E). (5.2) 

r.r' 

In eq. (5.1), the momentum of the center of mass co- 
ordinate (r -|- r')/2 is set to be 0, since we consider the 
spatial average. 

First, we consider the VIG case of Fig. 2(c) for 
rih = 1/16 = 0.0625 and 5 = |nh. Figures 11 shows 
the contour map of N{k,E) — ^^^jN^ik^E) along the 
symmetry lines in the original Brillouin zone. We also 
calculate N{k,E) for other nh cases, and obtain qual- 
itatively the same structure. While the energy disper- 
sion splits into A''-bands by the N-site periodicity, the 
spectral weight in 27r x 27r fe-space has the structure re- 
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Fig. 12. Spectral weight N(k,E) in 2n X 2n fc-space at fixed E in the VIC case of Fig. 11. E/t = -0.56 (a), -0.39 (b), -0.14 (c), 
(d), 2.02 (c), 2.27 (f). We see the valence band in (a)-(c), the mid-gap state in (c) and (d), the conduction band in (e) and (f). kx 
and ky are scaled by tt. U = S.dt, t' /t = —1/6, njj = 1/16 = 0.0625, S/n^^ = 3/4. The Fermi surface for <5/nij = 3/4 corresponds to 
the ridges in (d). The Fermi surface for (5/nh = 1/2 has the same structure as the ridges in (c). 



fleeting the original dispersion e{k) as shown in Fig. 11, 
since the amplitude of the wave function Uk^ ,j a,m for 
each m reflects the original dispersion e{k). Some mQ- 
shifted e(fc) by the band folding also appear in the figure. 
While their amplitude are usually small, they form the 
SDW gap and the mid-gap states. Darker region in Fig. 
11 means larger N{k, E). As the stripe is a ID structure, 
N{k, E) is two-fold symmetric. Then, there are differ- 
ences of N{k,E) between (1,0) and (0,1) directions in 
the VIC case. In Fig. 11, N{k,E = Ep) has little am- 
plitude around (1,0). But near (0, 1), the mid-gap state 
appears with large amplitude at Ep- There, we can see 
the split of the mid-gap state and the valence band be- 
low Ep. Around (^,^) and (— i,^), we see a gapped 
structure, where the mid-gap state does not appear at 
Ep. 

These features of the spectral weight are consistent 
with t|h£)se. -obtained by the ARPES experiments on 
LSCO.caOEa) Since there are two-type domains of the 
a;-direction and the y-direction stripe in the material, we 
observe the spectral weight as NY){k,E) — [N{k,E) + 
N(k^ ^ ky,E)]/2, i.e., the overlap of N{k,E) in Figs. 
11(a) and 11(b). 

When 5 is shifted keeping nh fixed, we confirm that 
N{k, E) is not qualitatively changed except for the posi- 
tion of Ep. With increasing i5, Ep is shifted to upper en- 
ergy. For larger 6 such as J = nh, the gap edge at (5,5) 
is located at deeper position below E-p. The ARPES 



experimentsEErtZl'tj-' reported that the gap at (5, ^) for 
a; = 0.15 disappears with increasing rih. This gap disap- 
pearance can be explained as follows jn. our scenario. In 
the elastic neutron scattering data,tlLl'lj'Ll' with increas- 
ing rih, 5 shows saturation and it is shifted from 5 ^ n\i 
io 5 '-^ ^n^. When 5 approaches i^h, Ep is effectively 
shifted toward lower energy, and touches the top of the 
valence band near {\:\)- Then, the gap in the (i, i)- 
direction disappears. As Ep approaches the gap-edge at 
(i, i) gradually, the gap disappearance is "observed" at 
larger 5 than 6 ~ ^^h, if the energy resolution is broad 
in the experiment. 

To clearly see the fc-dependent SDW gap, we consider 
N{k, E) in 2tt x 2ti fc-space at fixed E. It is shown 
in Fig. 12 for VIC. From (a) to (f), E is increased. 
In Figs. 12(a)-(c), the valence band appears. Its be- 
havior is understood by considering (fc) of the C- 
AF case in eq. (3.1). Below Ep, N{k,E) has ridges 
along the contour line for E = E-{k). Far from Ep, it 
reflects the original dispersion e(fc) [Fig. 12(a)]. For 
£;_(7r(l,0)) < E < £:_(7r(i,i)), N{k,E) disappears 
near (1, 0) and (0, 1), because the AF gap opens at these 
fc-points. With increasing E, the remained contour lines 
of ^^-(fc) decrease their length, and converge upon (5,5) 
[Figs. 12(b) and 12(c)]. Near Ep, we also see the frag- 
mentation of the ridges due to the x-direction ID stripe 
structure [Figs. 12(b)]. The valence band contribution 
disappears for E > i?_(7r(i, i)). 
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Fig. 13. The same as Fig. 11, but for the insulator DIG case of Fig. 1(b). U = 3.6t, t' /t = -1/6, = 1/16 = 0.0625, S/riy, = 1/2. 



The mid-gap state appears in Fig. 12(c) and 12(d). On 
raising E, it first appears near (0, 1) [and (1,0) with smaU 
amphtude] at the bottom of the mid-gap bands in addi- 
tion to the valence band contribution [Figs. 12(c)]. The 
"Fermi surface" in the VIC case corresponds to the ridges 
shown in Fig. 12(d), where the valence band contribution 
already disappears. Since the mid-gap state is locahzed 
along the a;-directional stripe, N{k, E) is distributed on 
Hues perpendicular to direction, like a ID Fermi sur- 
face. As it also reflects the original Fermi surface of e(fc), 
the amplitude of iV(fc, E) is large near (0, 1). Then in the 
spectral weight of the domain structure Nj){k,E), the 
Fermi surface exists near (0,1) and (1,0), but disappears 
near (i, i). This fe-dependent Fermi pirface structure is 
consistent with the ARPES results. EZfEa) We obtain qual- 
itatively the same Fermi surface also for 5 — n\y. It is 
qualitatively the same results of the case <' = under the 
test stripe potential in ref. 21, For 5 — ^rih, since E^ 



is decreased effectively, the Fermi surface has the same 
structure as shown in Fig. 12(c), i.e., the valence band 
contribution also appears near (^, ^) in addition to the 
mid-gap state's contribution near (1,0). This feature is 
contrasted with the 5 = |nh case as mentioned above. 

With increasing i?, the larger ridge-wall near (0, 1) 
shown in Fig. 12(d) goes away from the line kj. = 0, 
and the small ridge-wall near (1,0) approaches fc^ = 0. 
In this shift, the ridge of N{k, E) has large peaks near 
the fc-points of the original Fermi surface by e(fe). Then, 
if d is randomly distributed (it means effectively that 
Ef is randomly distributed), the original Fermi surface 
shape of e(fc) appears as an averageBI^ 

The conduction band contribution appears in Figs. 
12(e) and 12(f). It qualitatively follows E+{k) dispersion 
in eq. (3.1). For E+(tt{1,Q)) < E < £;+(7r(i, i)), the 
conduction band contribution appears only near (1,0) 
and (0,1) [Fig. 12(e)]. With increasing E, it spreads 
and approaches (^, ^), and it is connected each other. 
For E > -B+(7r(i, i)), the ridges of N(k, E) reflect the 
original dispersion e(fc). We also see the small amplitude 



of the circle ridge around (0,0), which is the Q-shift of 
the original e(fe) contribution. 

Next, we consider the DIG case calculated for the same 
parameter set. It is shown in Figs. 13 and 14. There, 
/Q-shifted e{k) dispersion curves form the valence band, 
the conduction band and the mid-gap state as in the VIC 
case. As shown in Fig. 13, a small gap opens between the 
valence band and the mid-gap state in the DlC case, and 
Ep is in the gap. Reflecting the diagonal-direction stripe, 
N{k, E) is two-fold symmetric. The differences appear 
between the (i, i)-direction and the (— i, i)-direction. 
In Fig. 13(c), N{k,E) shows flat dispersion near E-p 
along the line (—1,0) - (0, 1). It reflects the ID-like dis- 
persion due to the stripe, as discussed below. 

To clearly see the fc-dependent SDW gap for DIC, we 
consider N(k, E) at fixed E as shown in Fig. 14. The 
valence band contribution appears in Figs. 14(a)-(c). 
For £;_(7r(l,0)) <E < E^{tt{\, \)), N{k,E) near (1,0) 
and (0, 1) disappears as in the VIC case. Further, by 
the effect of the diagonal stripe, the remained ridge-wall 
near (i, i) is split into two peaks [Fig. 14(a)], and shifted 
toward (|,|) [Fig. 14(b)], and disappears [Fig. 14(c)], 
with increasing E. As for the ridge- wall near (— ^, ^), its 
length monotonically shrinks on raising E [Figs. 14(a) 
and 14(b)], but the ridge survives for a while after the 
(5, 5)-ridge disappears [Fig. 14(c)]. 

Figure 14(d) shows the mid-gap state's contribution. 
As the mid-gap state is localized along the diagonal 
stripe, its dispersion is ID-like. The ridge-wall ap- 
proaches the line k^ — ky with increasing E. The ridge- 
wall has two large peaks. Their peak-positions are near 
the contour line of the original dispersion e(fc). 

The conduction band contribution appears in Figs. 
14(e) and 14(f). At the bottom of the conduction band, 
the peaks of N{k,E) appear around (1,0) and (0,1) 
[Fig. 14(e)]. And the additional ridges connect these 
peaks along the line (1,0) - (0,-1). The peaks are sepa- 
rated along the (1,0)- (0, 1) direction. This ID-like ridge 
comes from the diagonal stripe structure. With increas- 




Fig. 14. Spectral weight N{k, E) in 2-k x 2-k fc-space at fixed E in the DIG case of Fig. 13. E/t = -0.56 (a), -0.31 (b), -0.14 (c), 0.27 
(d), 1.77 (e), 2.27 (f). We see the valence band in (a)-(c), the mid-gap state in (d), the conduction band in (e) and (f). and ky are 
scaled by tt. U = 3.6t, t' /t = -1/6, = 1/16 = 0.0625, 5/nh = 1/2. 



ing the contributions near (1, 0) and (0, 1) spread out, 
and they are connected each other. Then, N{k^ E) has 
a connected ridge along the contour hnc of the original 
dispersion e(fe) [Fig. 14(f)]. We also see the Q-shiftcd 
contribution of e(fc) as small ridges. 

§6. Optical Conductivity 

We calculate the optical conductivity by the inter- 
band vertical transition within the reduced Brillouin 
zone. Since there are N bands in the stripe structure 
with the iV-site periodicity as shown in Figs. 1 and 2, 
the inter-band absorption occurs by these N bands. The 
optical conductivity Recr^ ,j' {ijj)Ajfrd'r-^ y) is given by 
the Kubo-Greenwood formula.tSoo-' In our notation, 
it is written as 

RecTj- J' (uj) 

= ^ E /(^fco,a,.)[l-/(i^fco,a'..)] 

X (fco, a, a\Jj\kQ,a' , <j){ko, a, a\Jji \ko,a' ,cr) 

(6.1) 

Equation (6.1) is derived from the current-current corre- 
lation function. The x-component of the current opera- 
tor is defined as X{r,) = -i|e| 'EAti4+x,a ~ 4-x,a) + 

t'i4+x+y,a + 4+x-v,a - 4-x+v,a " 4-x-yM)}'^^.T in OUr 

square lattice case, where i + x and i + y mean {i^ + 1, iy) 



and {ix^iy + 1), respectively. Then, the a;-component of 
the oscillator strength in eq. (6.1) is given by 

{ko,a,a\Jx\kQ,a',a) = \e\^[t{uko^a,criri+x) 

i 

~'^ko,a,<ji.^i—i)} ^ {'^ka,a,ry{xi+x+y) + '^kf),a,<ji.^i+x — y) 
~'^ko,a,<j{''^i—x+y) ~ ^fco ,a,c7 (''■i-S-y ) }] ^fco 

(6.2) 

instead of the usual form (V-'fcl^xlV'fc')- 

Figure 15 shows Re(T||((jj) for the parallel field to the 
stripe and Kea^luj) for the perpendicular field in the 
VIC case. Here, we do not consider the Drude absorp- 
tion by the intra-band inelastic scattering. For rih — 
1/26 ~ 0.0385 [Fig. 15(a)], there is a large peak near 
Lo ^ 2.5t. It comes from the valence band conduction 
band absorption with larger to than the AF gap, and 
survives also in the C-AF state. We call the peak as the 
"AF peak" . The peaks a,t uj/t ~ 1 ^ 1.5 comes from the 
absorption related to the mid-gap states. We call them 
as the "mid-gap peaks" . They vanish in the C-AF state. 
For larger rih [Fig- 15(b)], the AF peak is suppressed and 
the mid-gap peaks are enhanced. The mid-gap peaks of 
Kea±{Lu) are slightly shifted to lower ui compared with 
RetT|| (uj). 

The Die case is shown in Fig. 16. With increasing 
nh, the AF peak is suppressed and the mid-gap peaks are 
enhanced, as in the VIC case. There are large differences 
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Fig. 15. Optical conductively Re(T|| (oj) [solid lines] and Hca±(w) 
[dashed lines] for VIC in an arbitrary unit. Recu (u;) (Reaj_{Lu)) 
is for the field parallel (perpendicular) to the stripe. The Drude 
absorption is not included in this figure, (a) nh = 1/26 = 0.0385, 
5/nh = 13/18. (b) rth = 1/8 = 0.125, (5/nh = 5/7. U = 3.6t, 
t'/t = -1/6. 



between RetT||(ijj) and Ref7j^(cj) in the DIG case. The 
mid-gap peaks of Recr^ (w) are large, and they are located 
at larger uj. For rth = 1/8 = 0.125, the AF peak vanishes 
for Ke(j±{uj), but still has large height for Re(7||((jj). 

The experiments of the optical conductivity were per- 
formed on LSCO [ref. |§ and La2-a;Sr^Ni04 [ref. |l)- 
There, the broad peak develops in the mid-infrared re- 
gion with increasing nh- It may be related to the mid-gap 
peak's contribution of the stripe state. 

§7. Summary and Discussions 

The stripe structure in high-Tc cuprates and its elec- 
tronic structure are studied by the self-consistent mean- 
field theory of the Hubbard model. The SDW gapped 
insulator is changed to a metal by introducing the re- 
alistic Fermi surface topology in VIC. For appropriate 
U, DIG is changed to VIC at critical rih with increasing 
rih- The DIC is insulator, and VIC can be metallic. In 
the insulator, 5 = ^nh. But, in the metallic VIC, S can 
be larger than in^. These features seem to be consis- 
tent with the exj^imental results of the elastic neutron 
scattering,l§ByBBUEi) as discussed in §| 

The SDW and CDW of the stripe structure are a si- 
nusoidal wave at large nh- They become, respectively, 
square and solitonic waves at small nh- There, excess 
carriers are localized in the stripe region and the outer re- 
gion is half-filled AF. The localized state forms the mid- 
gap state inside the AF gap. We calculate the LDOS, the 
spectral weight and the optical conductivity, and clar- 
ify the contribution of the mid-gap state. In the metal- 
lic VIC case, only the mid-gap state appears at Ep for 



Fig. 16. Optical conductively Re(T||(t<j) [solid lines] and Ilca±{uj) 
[dashed lines] for DIC in an arbitrary unit, (a) njj = 1/26 = 
0.0385. (b) nh = 1/8 = 0.125. U = 3.6i, t'/t = -1/6, <5/nh = 
1/2. 

6/n\i ^ I or larger. Then, while the stripe region is 
metallic locally in space, the outer AF region is an insu- 
lator. IVIetallicity attains locally on the stripe region for 
finite rih- In the spectral weight, the contribution of the 
mid-gap state has large amplitude near (0, 1) and little 
amplitude near 

(1,1) at Ep in VIC- Then, the "Fermi 
surface" disappears near ^), leaving the Fermi surface 
"arcs"- On the other hand, for S/uh ~ ^, both the mid- 
gap state and the valence band appear at Ep. In this 
case, both the inside and outside spatial regions of the 
stripe are metallic- In the spectral weight, the valence 
band contribution appears near (5,5) in addition to the 
mid-gap state's contribution near (0, 1). Then, the Fermi 
surface is observed as a connected one. These features 
of the spectral-weight are consistent with the ARPES 
experiments,^'^^^^ as discussed in i|5[ As a result, the 
electronic structure near Ep is determined by the ratio 
of the incommensurability and the hole density, S/rih, 
which is a key parameter of the problem. We expect 
that the features of the electronic state obtained here 
are examined by the experiments such as ARPES, STIVE 
and the optical conductivity. These experiments may be 
powerful methods to obtain vital information about the 
mechanism of the stripe structure in high-Tc cuprates. 

So far, the metallic stripe state with 6/nh 1 could 
not be reproduced by the mean field theory, which gives 
an insulator with (5/nh = 5- Then, the stripe state was 
said to be a anomalous state where the mean field theory 
can not be applied. However, by introducing the realis- 
tic Fermi surface topology, the mean field theory can 
gives the metallic state with S/n\i = 0.75 ^ 0.88 (which 
approaches 1 depending on the Fermi surface topology) 
even in a simple Hubbard model. While the mean field 



iviasanori ichioka ana ivazusnigc iviachiua 



theory over-estimates the stable parameter region of the 
ordered states or the AF moments, the obtained features 
are quaUtatively vaUd as a first approximation within the 
ordered state. Then, the mean-field theory can be a first 
step for further studies to consider the effects of the extra 
interactions or fluctuations, and to obtain the behavior 
of the various physical quantities in the stripe state. 
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